Impurity systems fluctuating between two magnetic configurations 
candidates for non Fermi-liquid behavior 
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The appropriate generalization of the isotropic impurity Anderson model for va- 
lence fluctuations between two magnetic multiplets l n and l n+l is solved in the strong- 
coupling limit of Wilson's renormalization group for I < 3. Except in the extreme 
case of j — j coupling, the ground state is degenerate, the impurity magnetic moment 
is very small and is overscreened by the conduction electrons. The strong-coupling 
fixed point turns out to be unstable. Thus, at low enough temperatures the physics 
of the systems should be governed by an intermediate coupling fixed point and non 
Fermi-liquid behavior is expected. 
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The Anderson model for dilute magnetic alloys and its integer valent limit, the Kondo 
problem, have played a decisive role in the understanding of intermediate- valence and heavy- 
fermion systems. The scaling behavior of these models has been established [|l]-|4j] . The low 
temperature properties are described by the strong-coupling fixed point, characterized by a 
non-magnetic ground state and Fermi-liquid behavior. These models and their extension to 
an arbitrary degeneracy of the magnetic configuration have been solved exactly using the 
Bethe ansatz ||. The predicted thermodynamic properties are in general in good agreement 
with experiments in dilute systems fluctuating between one magnetic and one non-magnetic 
configuration ||. 

The properties of systems in which both accessible configurations are magnetic are 
markedly peculiar. For example TmSe is the only intermediate valence compound which 
orders magnetically with a local magnetic moment 1.5/ie 0, much smaller than those of 
the 4f 12 and 4f 13 configurations. The compounds UPt3 and UBei3 || also order magneti- 
cally with a magnetic moment of the order of 10~ 2 /i£. The magnetic susceptibility of dilute 
Tm systems points towards a magnetic ground state of Tm impurities |5|]. The magnetic 
susceptibility of the Y0.gU0.2Pd3 also diverges at temperature T — > and a finite entropy at 
T = 0, S(0), is suggested by specific heat measurements |r0yiT|| . D.L. Cox stated that the 



different behavior of U and Ce systems might be caused by different underlying symmetries, 



and constructed an Anderson-like model starting from crystal- field split levels |E|]. In the 
limit in which only the T 3 non-magnetic doublet of the 4^ configuration and the T 7 states 
of the excited 4f* configuration are retained and also the conduction j = 7/2 partial- wave 
states are omitted, his model takes the form of the spin 1/2 two-channel Kondo problem 
13| . The low-temperature properties of the overscreened multi-channel Kondo problem are 



characterized by an intermediate-coupling fixed point [13, Tm with a no n- Fermi-liquid be- 



havior [|14|,[nj which might be related to the marginal Fermi liquid phenomenology of the 



high Tc systems [|TIJ. The spin 1/2 two-channel Kondo problem has been solved exactly (15 
and is consistent with several properties observed in Y0.gU0.2Pd3, such as the logarithmic 
divergence of the specific heat at low temperatures with a residual entropy 5*(0) ~ (1/2) ln2 



T0|| - However, other properties, such as the magnetic field dependence of the specific heat, 



are inconsistent with the model |L1]] . The same happens with some data of UBei3 [[TIJ . 

The study of microscopic models with non-Fermi-liquid behavior, particularly if they 
are relevant to real systems, is of great interest at present. In this letter we show that 
neglecting crystal fields, the low-temperature properties of systems fluctuating between two 
p, d or f configurations are characterized by an intermediate-coupling fixed point, except in 
the unrealistic extreme cases in which one of both hybridization channels Vj(j = I ±1/2, see 
below) vanishes. Unfortunately, only one Vj was retained in most previous studies of valence 
fluctuations between two magnetic configurations |T7| - |2TJ| . In particular V/ + i/2V$_i/2 = 0, if 
j — j coupling (instead of the more realistic LS coupling) is assumed. Also, as we have 



shown elsewhere ||22|| , there are not enough conduction electron degrees of freedom in the 



variational wave functions considered in Ref. |2TJ, to capture the essential physics of Tm 
impurities. 

The generalized Anderson model for an impurity ion in an isotropic medium can be 
written in the form: 

H ttband H m i x -\- Hion (1) 

The first term describes a band of extended states and we need to retain only states with 
the same orbital angular momentum as that of the impurity shell / 

Hband — e k C \ncj C kfjL(T = X/ e k C kjm C kjm (2) 
kfia kjm 

where /i, a are the orbital angular momentum projection and spin respectively. The last 
expression is an alternative description in terms of the total angular momentum j = I ± 1/2 
and its projection m. Similarly: 

H mix = Vj2(flc ka + h.c) = V^2{f] m c kjm + h.c.) (3) 

kfia kjm 

where the / operators refer to impurity electrons. Hi on describes the impurity shell and 
contains all the correlations of the problem. We retain only the ground state multiplets of 
two neighboring l n and l n+1 configurations: 



H lon = EJ2 I JoM >< J M | +{E + A) I Ji^i >< ^1^1 I (4) 

M Mi 

where Jj, Mj are the total angular momentum and projection of the l n+t configuration. 
We find more convenient to express the / operators in terms of the Hubbard operators 
| JqMq >< J 1 M 1 |. When this is done H mix takes the form: 

H mix = V i< J ojM m | J X M X > (4 im | J M >< J X M X \ +h.c.) (5) 

kjmMoMi 

where the < JojMom \ J\M\ > are Clebsch-Gordan coefficients. The Vj/V are not free 
parameters, but are determined by the particular form of the highly correlated ionic states. 
In the simplest case of j — j coupling, it is easy to verify that for < n < 21 — 1: 

V l+1/2 = 0; = V< J {1 - 1/2) J {1 + 1/2 - n) \ J X J X >^ (6) 

while for 21 < n < Al + 2: 

VU/2 = 0; V l+1/2 = V < J (l + 1/2) J (3Z + 1/2 - n) \ M >~ x (7) 



Instead, using LS coupling and the Hund rules, it can be shown that 

Vj _ E M< ri?s? < L S L Z S^ | J J >< LxSiLlSl \ M >< l\[io \ j(J x - J ) >< L IL z /j \ L X IL\ >< Sq^S, 
V < J jJo(Ji - Jo) | JiJi >< L lL (Li - L ) I LiLi >< S , |5 , (5'i - S ) \ S 1 S 1 > 

(8) 

The result of Eq. (8) for the cases in which both l n and l n+1 configurations are magnetic is 
given in Table I. 

The integrability of the model using the Bethe ansatz has been studied by Aligia et al. 



| 17| . When both configurations are magnetic, only the cases with / = 1 (see Table I) are 



integrable if in addition V3/2 = is assumed. They were solved by the same authors ]I7[ and 



by Schlottmann [[J. In the integer- valent limits | A |>> V, these cases are equivalent to 



underscreened Kondo problems (T7j , characterized by the same strong coupling fixed point at 
low temperatures as the original model . The ground state has total angular momentum 
min(Jo, Ji)- It has been shown recently that in the underscreened Kondo problem, Fermi- 
liquid excitations coexist with the magnetic degrees of freedom and that they are decoupled 
at low frequency |[25||. 



Unfortunately, for the most realistic cases for valence fluctuations between two magnetic 
configurations the model is not integrable with the Bethe ansatz and also alternative, accu- 
rate enough methods used for I = 0, have serious practical limitations. The nature of the 
correlations of Hi on would require several different Hubbard- Stratonovich transformations 
for the implementation of a Monte Carlo algorithm [EH], while in Wilson's renormalization 



group []2|-§J , the size of the relevant Hilbert space is enlarged by a factor 2 4/+2 in each in- 
teration, severely affecting the numerical accuracy for / 7^ 0. A similar, although not so 
severe restriction was encountered in the renormalization group treatment of the spin 1/2 



two-channel Kondo problem [[14]] . However, formally, Wilson's renormalization group can be 
implemented in the same way as for I = |2|-|J. The only difference is that the logarithmic 
discretized band has orbital degeneracy. We call c n ^ a the states and destruction operators 
of the different onion-like spherical Wannier functions (f nfl in the notation of Ref. 3). It is 
clear that if one starts from a small value of V, during the first few iterations the system 
is near the free-orbital fixed point and the effective value of V increases. If in the subse- 
quent renormalization group flow the system does not reach a fixed point at an intermediate 
value of V, then, the low temperature behavior is described by the strong-coupling fixed 
point V — > 00, as in the case of / = PH|] or a particular model for valence fluctuations 



between two magnetic configurations p4| . However, if the strong-coupling fixed point is 
unstable, there must exist at least one intermediate-coupling fixed point, as in the spin 1/2 
two-channel Kondo model fl3 , |14" l. In this case, there are no asymptotic scattering states 



14| and a non-Fermi-liquid behavior is expected. 

Thus, the study of the strong-coupling fixed point might provide a clue about what 
kind of physics is expected in the model at low temperatures. The properties of this fixed- 
point are determined by the local Hamiltonian H in which only the innermost Wannier 
functions c 0fl(T are considered and all energies are neglected compared to V [EJ-[| . This one- 
site problem is far from trivial for realistic systems. The size of the Hilbert space of Ho is 
(Jo + J\ + 1)2 4?+3 . We have found the ground state of H in each subspace of a given number 
of particles n + n c and total angular momentum projection J*, using the Lanczos method 



and a modified version of it ||27|| . The quantum numbers of the ground state are given in 
Table I. They are independent of A and Vi±i/ 2 as long as both Vj 7^ 0. The cases with 
the same I and Jo and J\ interchanged are related by a special electron-hole transformation 
fi7| , |22|j . In general, the ground state is highly degenerate. There are states differing not 
only in J*, but also in J t and n c with the ground state energy. Exceptions are the cases 
with I = 3 and 1 < n < 4 where only one value of n c and Jt characterizes the ground state. 
However, due to numerical uncertainties, we cannot completely rule out the possibility that 
states with other values of Jt are degenerate with the ground state in these cases. 

The expectation value at the impurity magnetization J? = Si,Mi M I JiMi > is always 
negative for any state of the ground state with maximum J t z . Some specific values are 



given in Ref. |[28|| . Thus, the impurity magnetization is overscreened. The absolute value of 
< Jf mv > decreases when the magnitude of the Vj are changed from the ideal LS'-coupling 
values to those corresponding to intermediate coupling, and vanishes in the j — j coupling 
limit. In this limit the degeneracy of the ground state becomes irrelevant, since the states 
c jm with j such that Vj = become decoupled. If one is almost in the j — j coupling limit, 
say Vj\ » Vj2, where Vj2 would vanish according to Eqs. (6), (7), there will exist a regime 
at intermediate temperatures characterized by a fixed point in which Vj\ = oo, Vj2 = 0. 
Our results for this fixed point agree completely with those obtained by Lustfeld for / = 3 
and several values of n in the integer valent limit of the present model 0]. In particular 
the ground state is non-magnetic, in agreement also with Refs. [ I9p0 . However, this fixed 
point is unstable if Vj2 > 0, and a finite though small impurity magnetic moment as found 
in T m and U |7]|| compounds starts to develop at lower temperatures. 

For 1=1 and V3/2 = we find a ground state with J t = min(Jo, J±) in agreement with 
Bethe ansatz results and renormalizat ion-group calculations . 

In the following we discuss the stability of the strong-coupling fixed point. In the imme- 
diate neighborhood of this fixed point, eliminating the potential scattering by means of the 
mapping explained in Ref. P] Appendix A, section 4, the Hamiltonian has a form like: 
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H = H (V) + H' band + YA*°*> + h - c -) ( 9 ) 

where H (V) is the problem solved above with a very large but finite V and H' band represents 
the band without the cq^ states, discussed in Section III.B.3 of Ref. 0. In most of the 
present cases in which different values of n c are present in the ground state, the operators 
Co^o. and Co M(7 connect different states with the ground state energy. In each renormalization 
group iteration there is, by definition, an overall factor A 1 / 2 where A > 1 is the parameter 
of the logarithmic discretization of the band (Eq. (2.22) of Ref. 2). Also the operator 
introduces a factor A 1//4 in each iteration (Eq. (4.23) of Ref. 2). Thus, the last term of 
Eq. (9) is accompanied by a factor A 1 / 4 and is a relevant operator, indicating that the 
strong-coupling fixed point is unstable. 

The cases in which only one n c is present in the ground state (1=3, n=l or 2 in Table 
I) are more delicate, but are completely analogous to the corresponding case of the spin- 
1/2 two-channel Kondo model |[3[|14| . In these cases, the last term of Eq. (9) affects the 
eigenvalues of H in second order perturbation theory, leading to an effective interaction H e f f 
between the ground state of H and the Ci^ a states. Following a similar reasoning to that 
which led to the effective exchange Hamiltonians in the integer valent limits of the present 
model [17,1^], one can show that H e ff has the form of an exchange Hamiltonian between 
the total angular momentum J t of the ground state of Hq and the ci jm states. Also since the 
relevant denominators in the perturbative treatment should increase like V as A 1//2 on each 
iteration, H e ff is a marginal operator. Following Ref. 2, we expect that the strong-coupling 
fixed point is stable if the coupling constant J e ff of H e ff is negative (ferromagnetic) and 
unstable otherwise. A simple argument due to Nozieres and Blandin shows that if the 
impurity magnetic moment is over-(under) screened in H , J e ff > (J e ff < 0). Thus, the 
strong-coupling fixed point is also unstable. 

In conclusion, valence fluctuations between two magnetic l n and l n+1 configurations, 
except in limiting unrealistic cases in which Vj+i/2T^-i/2 = (e.g. extreme j — j coupling) 
are epected to be characterized by an intermediate coupling fixed point at low temperatures, 
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in absence of crystal field. Although the cases with / = 1 might look academic, they are the 
simplest ones to be studied by alternative methods, which we hope will follow the present 
study. The unrealistic case I = has already provided a great deal of understanding of the 
physics of systems fluctuating between one magnetic and one non-magnetic configurations. 

We would like to thank Andrei E. Ruckenstein for many enlightening discussions, and 
L.M. Falicov for a critical reading of the manuscript. One of us (MB) is supported by the 
Consejo Nacional de Investigaciones Cientificas y Tecnicas (CONICET). AAA is partially 
supported by the CONICET. 
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TABLES 



Parameters and ground-state quantum numbers in the strong coupling limit, for valence 
fluctuations between any two magnetic Hund-rules ground-state multiplets of neighboring 
l n and l n+l configurations with I < 3. Jj is the total angular momentum of the l n+l con- 
figuration and Vj is the hybridization energy for conduction electrons with total angular 
momentum j, using LS coupling. The numbers in brackets are those Vj which vanish for 
j — j coupling. The total number of particles and the total angular momentum of the local 
Hamiltonian H , are denoted by n + n c and J t respectively. Some cases of marginal interest 
with I = 3 were not solved. 
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(0.8165) 
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